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Abstract 

Based on a considerable amount of experimental evidence for generic properties of lateral organi- 
zation of lipid membranes containing inclusions, we first present a general model system of bilayer 
membranes embedded by nanosized inclusions, and account well for a series of unexpected behaviors 
in related experimental findings. (1) The appearance and disappearance of lipid/inclusion-rich rafts 
are observed with increasing the inclusion content. (2) The chain arrays of inclusions may form at 
high concentrations. (3) Location of inclusions changes with increasing the inclusion content, and 
may undergo a layering transition from one-layer located in the center of the bilayer to two-layer 
structure arranged in opposing leaflets of a bilayer. (4) The membrane fluidity is enhanced by the 
presence of inclusions. Our theoretical predictions address the complex interactions between mem- 
branes and inclusions, suggesting a unifying mechanism which reflects the competition between the 
conformational entropy of lipids favoring the formation of lipid-rich rafts and the steric repulsion 
of inclusions leading to the uniform dispersion. The present study advances our understanding 
of membrane organization by unifying these experimental evidences of real biomembranes with 
inclusions which can be different, but with the hydrophobic and rigid properties. 
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Recently, there is growing evidence that due to the presence of inclusions within mem- 
branes, the distribution of lipids is inhomogeneous, where lateral segregation could induce 
the formation of lipid/inclusion-rich raft domains. For instance, cholesterol which is one of 
the most important regulators of lipid organization, prefers to have conformationally ordered 
lipid chains next to it due to its hydrophobically smooth and stiff steroid ring structure 
and promotes the formation of lipid rafts(see, for exaniple, reviews[l|, ^, y] and recent re- 
search works^). On the other hand, it was reported 6] that hydrophobic drugs such as 



taxol (paclitaxel) 



and dipyridamole (DIP) [8]) may assist the formation of lipid/drug- 



enriched raft domains, and increase with increasing the drug content, but disappear at high 



concentrations, which has also been observed in cholesterol-lipid systems 



the perturbed lipids may lead to chaining of cholesterols 



or drugs i 



9|. Particularly, 



111 ] inside the 



bilayer. Furthermore, the introducing of taxol drug into the lipid layer will perturb the 



hydrocarbon chain conformation, which may increase membrane fluidity 0, 3, III- Most 
recently, some foreign inclusions such as silver nanoparticles were reported to have similar 
effects on the membrane fluidity 12]. 

Despite the common properties of membrane organization due to distinct inclusions and 
despite its important significance in cellular functions such as signal transduction and mem- 
brane trafficking [3|, the influence of the inclusions on such a change in lateral organization 
has not yet been considered in computational and theoretical investigations, and further 
insight into the mechanisms behind general evidence from lipid-inclusion complexes remains 
poor [13]. Previous theoretical works were concerned with the hydrophobic mismatch interac- 

15 . 16 ] and the possible formation of lipid rafts[l|, [3]. However, 



tion between inclusions 
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a detailed structural change with varying the inclusion content has not been systematically 
elucidated. In this letter, we examine a simple model that not only allows us to present 
a unifying description of these phenomena with varying inclusion concentrations, but also 
sheds light on physical mechanism behind membrane organization due to distinct inclusions 
such as intrinsic membrane protein, cholesterol, hydrophobic drug, or other bionanoparticles. 

Consider a lipid bilayer membrane containing Ud inclusions of radius R in an aqueous 
environment (Fig. 1). The volume of system is V = L^. X Ly X Lz, where and Ly 
are lateral membrane lengths under periodical boundary conditions, and is the size of 
system along the membrane normal direction. The bilayer membrane is composed of one 
type of lipids with two hydrophobic tails. The number of lipids is given by = 2 x 
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FIG. 1: (color) A schematic of the hpid bilayer membrane containing inclusions. 
a X Lj. X Ly, where a is the number density of hpids in one leaflet. The headgroup of 
lipids has the volume Vh, and two equal-length tails composed of N/2 segments each are 



assumed as flexible Gaussian chains 



17( 1 . The segment has the volume ^ and the length 



a. Thanks to the smooth and rigid property of bioinclusions compared with the flexible 
lipids , we assume inclusions as hard-sphere nanoparticles. The concentrations of 

the head and tails of lipids are (ph = niVh/V and (pt = niNpQ^ /V , and the inclusions 
(pd = AirR^rifi/SV. The solvent molecules have the volume Vs and the concentration 0^ = 



1 



18| . Recently, the self-consistent field theory (SOFT) has 



jeen proven 



to be powerful in calculating equilibrium morphologies in polymeric systems 



19 
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2l|, 



22 1 to account for 



while nanoparticles can be treated by density-functional theory (DFT) 
steric packing effects of particles. Interestingly, Thompson et al. 191] developed a hybrid 
SCFT/DFT approach to study mixtures of diblock copolymer and nanoparticles. On the 
other hand, the SCFT method is extended to successfully study the phase behavior of pure 



lipid systems 



13 
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231]. Here, we will extend the hybrid SCFT/DFT approach 191] to 



calculate the structural organization of the bilayer membrane in the presence of inclusions. 
The resulting free energy F for the present system 2J] is given by 



NF 



PoksTV 



M 



Qi 
V<Pt 



^ln(^ 

Ots V(i)s 



0d j^^QdO, 

a Vcpd' 



+ ^ j dr[xthNift{r)ifh{r) + XtsNift{r)Vs{r) 
+XtdNipt{r:)ifdijc) +XhsNiphir)ipsijc) 
+XhdNiph{Tc)^d{Tc) +XsdNips{r)ipd{r) 
-Wt{r)ipt{r) - Wh{r)(ph{r) - Wsir)(ps{r) 
-Wd{r)pd{r) - ^(r)(l - i^t{r) - i^h{r) 



(1) 



where a, = Vspo/N, and a = AirR^po/^N. Xth, Xts, Xtd, Xhs, Xhd, and Xsd are the Flory 
interaction parameters between tail-head, tail-solvent, tail- inclusion, head-solvent, head- 
inclusion, and solvent-inclusion, respectively. 9?t(r), v^hli"), V^(i(r), and (ps{r) are the local 



FIG. 2: (color) The pure bilayer membrane in an aqueous environment, (a) The top view of 
averaged density distribution of hpid tails in the upper leaflet, (b) laterally averaged density 
profiles of solvent((/?s), lipid tail{ipt), and head group((/j/i) along z axes, (c) The x-z cross-section of 
density distribution of lipid heads averaged in the y-direction. 

volume fractions of lipid tail, head group, inclusion, and solvent, and Wt(r), Wh{r), Wd{r), 
and Ws{r) are corresponding self-consistent fields, .^(r) ensures the incompressibility of the 
system, and Pd(r) stands for the inclusion center distribution. Qi, Qg, and Qd are respective 



191 ]. The last term in Eq. 



1) is the 
m. The 



partition functions for lipid, solvent, and inclusions 
nonideal steric interaction term 



22] with the weighted inclusion density ^^(r 
fields and densities are then determined by minimizing the free energy in Eq.(l), and the 
resulting self-consistent equations can be solved numerically 191]. To reasonably describe the 
interactions of hydrophobic inclusions dispersed in bilayer lipids with hydrophilic heads and 
hydrophobic tails, we choose XhsN = 0, XtdN = 3.5, XhdN = 15.0, XtsN = 15.0, XsdN = 
20.0, and XthN = 25.0. The other parameters are fixed to be = 30, Vh = Vs = Gpo^, 
a = 0.08, Lx = Ly = 60a, and = 15a. Here, the chosen was large enough, ensuring 
the solvent concentration 0^ = 1 at 2; = and L^. 

We first examine the case of a bilayer membrane in an aqueous environment in the absence 
of inclusions. Figure 2a shows the lateral distribution of lipid tails in one leaflet, which is 
uniform. The distribution is symmetrical with respect to the opposing leaflet. Figure 2b 
provides the average distribution of headgroups in the x-z cross-sections, reflecting that the 
shape of membrane surfaces is smooth in the absence of inclusions, since the lipid length is 
matched. In this case, the membrane is in a gel phase for the limited lateral mobility |9l. Ill|. 
where lipid tails tightly pack and hardly move in the pure membrane. Figure 2c shows 
the laterally averaged density profiles of ipt, fh, and ips across a bilayer, which displays the 
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FIG. 3: (color) The top views of averaged density distributions of inclusions (left) and lipid tails 
(right) in upper leaflet, (a) (/>rf = 0.10, (b) = 0.15, (c) (/>d = 0.20, (d) (/>d = 0.24, and (e) 

= 0.32. Here we fix = 2.5a which is smaller than the membrane thickness, 
basic structure of the bilayer and is favorably comparable with mesoscopic simulation and 
other coarse-grained models describing the bilayer composition 23|, |25|]. This validates the 
used SCFT which can reasonably explore conformational properties of lipids, in contrast 
to phenomenological models that ignore much of the internal structure of the bilayer 13| . 



Particularly, the approach will become powerful in exploring the lateral inhomogeneity of 
membrane with inclusions, in contrast to other field theory or simulation schemes[l3[], where 
the only one-dimensional density profiles along the membrane normal is shown in the cost 
of including molecular details of lipids. 

Figure 3 shows the in-plane density distributions of inclusions in the left column and the 
lipid tails in the right column with increasing the inclusion concentration 0^- When (l)d is low. 
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FIG. 4: (a) The translational entropy Sd of inclusions, (b) The translational entropy Str of hpids, 
(c) the conformational entropy Sc of lipids, and (d) the steric energy Fgt of inclusions vs the 
inclusion concentration. 

the translational entropy of inclusions has a significant contribution to the free energy of the 
system. Any compositional fluctuation of inclusions may lead to the lateral inhomogeneity 
of lipid composition, in which lipids are depleted in the inclusion-rich regions for ensuring 
the large translational entropy of inclusions(Fig. 3a). Astonishingly, as (pa is increased, 
the lipid/inclusion-rich rafts appear(Fig. 3b). In this case, the entropic contribution of 
lipids becomes significant, and tail chains clo sing to the rigid surfaces of inclusions can 



get extra conformational flexibility Q, ll6|, 126|, which enriches the lipids surrounding 



inclusions. Therefore, inclusions are accumulated in certain membrane regions, which leads 
to the formation of lipid-rich rafts. The raft size may be enlarged by increasing (pa (Fig- 
3c). Interestingly, as 0^ is added to a certain range(0rf = 0.24 ~ 0.28), lipid-rich rafts 
disappear, but instead the uniform distribution of both inclusion and lipid tail occurs (Fig. 
3d). This unexpected behavior provides a strong support for the experimental findings 
in drug- membrane js] and cholesterol-membrane complexes. This is due to the strong 
steric repulsion from large numbers of inclusions, which leads to the uniform dispersion of 
inclusions. Further increase of (pd leads to the deformation of lipids which produces the 
effective attraction between inclusions. The deformed conformational entropy can partially 
be released by chaining of inclusions under the lipid-mediated attraction, shown in Fig. 3e. 
Such a regularly modulated inclusion-rich stripe structure has been reported in the drug- 
membrane gI, 11] and cholesterol- membrane complexes j^, [lO| where the lipids form ribbons 
between the aligned cholesterol domains. 

Figure 4a shows the translational entropy Sd of inclusions. For low (pd, the large increase 
of Sd with (pd indicates that Sd plays an important role at first, which can account well 
for the formation of the weak inhomogeneous distribution of lipids(Fig. 3a) to ensure large 
translational entropy of inclusions. For large (pd, Sd decreases with the appearance of the 
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FIG. 5: (color) Laterally averaged density profiles of solvent ((/^s), lipid tail((^f), head group(99/i), 
and inclusion (y?^) along z axes (left), and the x-z cross-sections of density distributions of lipid 
heads averaged in the y-direction(right). Parameters are the same in Fig. 2. 
chaining structure. Figure 4b shows the translational entropy Str of lipids. The curve goes 
up monotonously, meaning that the addition of inclusions increases the lateral membrane 
fluidity, which is a characteristic of gel-liquid transition [ol, The reason is that lipid 

tails are stretched along the bilayer normal, which increases translational degrees of freedom 
and thus enhances the lateral mobility of lipids Figure 4c shows the conformational 
entropy Sc of lipids. Beginning from a low conformational entropy in the gel phase of the 
pure membrane, lipids get more conformations with the perturbation of added inclusions 
where the lipid-enriched rafts are formed. However, there is a small decrease in the range 
(pd = 0.24 ~ 0.28 with the disappearance of lipid rafts, where all lipids are strongly stretched 
with the same length. Figure 4d shows the steric repulsion energy Fgt of inclusions, which 
increases with increasing 0^. Therefore, it is the steric repulsion between inclusions that 
suppresses lipid rafts. For high (j)d, the repulsion of inclusions is stronger, but at the same 
time, the deformed lipid-mediated attraction between inclusions becomes also stronger. The 
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only way that the conformational entropy of deformed lipids is released(Fig. 4c), is to drive 
inclusions to assemble anisotropically along one direction. As a result, the deformed lipids 
provide an additional lateral anisotropic interaction between inclusions, which stabilizes the 
parallel chain arrays of inclusions. 

Finally, Fig. 5 shows laterally averaged density profiles of different components across the 
bilayer in the left column and average density distributions of lipid heads in x-z cross-sections 
in the right column. For low (f)d{Fig. 5a), the inclusions assemble in the bilayer midplane for 
the membrane stability, which is displayed by one peak of density profiles {(fa) of inclusions. 
With increasing 0^, one peak density profile may be saturated (Fig. 5b), and further increase 
of inclusions leads to the occurrence of two peaks of (fd (Fig.5c-e), implying the relocation 
of inclusions where the two-layer distribution of inclusions is arranged in opposing leaflets 



of a bilayer 



261 ]. Previous experiment 



111 ] on drug-membrane complexes has also shown that 



the location of drugs in the bilayer depends on drug content. The two-layer distribution of 
inclusions ensures that the ends of lipid tails can still remain in the membrane midplane 
favoring the conformation of lipids. On the other hand, in Fig. 5a-c, the irregular density 
distribution of headgroups originates from the lipid chain-length mismatch, while in Fig. 
5d, the membrane surfaces become flat with a matched length of strongly stretched lipids 
where lipid-rich domains disappear. By comparison of the distances between two peaks of 
head profiles (iph) in the left side of Fig. 5, we also find that the thickness of membrane 
continuously increases with the addition of inclusions. 

This work was supported by the National Natural Science Foundation of China, Nos. 
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